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Abstract 

We construct a generic mass model of a disk galaxy by combining the usual two 
mass components, an exponential thin stellar disk and a dark halo for which we 
alternatively assume one of the two most common mass distributions. We construct 
its generic rotation curve, with two free parameters linked to its shape and one to 
its amplitude, namely the halo characteristic scale lenght, the fractional amount of 
Dark Matter and the total mass at one disk scale lenght. 

We find that there is no need of a fine-tuning in the above parameters to produce 
smooth, featureless rotation curves like the ones actually observed. The rotation 
curves obtained from different mass distributions are significanlty affected by the 
variation of the parameters, and show substantial difference in their slope; however, 
they do not feature irregularities marking the transition from a disk-dominated to 
a halo-dominated regime. 

On the other hand we show, by direct analysis of a large number of plausible and 
different mass models, that an observed rotation curve fulfilling well defined quality 
requirements can be uniquely and properly decomposed in its dark and luminous 
components; through a suitable mass modelling together with a maximum-likelihood 
fitting method, in fact, the disentangling of the mass components is accomplished 
with a very high resolution. 
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1 Introduction 



The rotation curves of spiral galaxies are observed to be smooth, featureless, 
out to several disk scale lenghts: they are essentially linear functions of the 
galactocentric distance, with a hint of flattening at the outermost radii (Per- 
sic, Salucci & Stel, 1996; see Appendix A). Rotation curves are determined by 
the full distribution of the mass of a galaxy, including its luminous component 
(the disk) and its dark component (the halo). While the radial distribution 
of the luminous matter is observed for spiral disks, the disk mass itself is un- 
known. Regarding the Dark Matter, both the amount and the distribution 
are unknown, and since these physical quantities affect the rotation velocity 
profile in a different way at different radii, the study of the rotation curves 
becomes a very important tool to investigate the structure of galaxies. 
There is sometimes the idea that the smoothness of the observed rotation 
curves implies a fine-tuning of the halo/disk structural parameters. Related 
to it, although logically orthogonal, there is the view that, due to this smooth- 
ness, it is impossible to disentangle the dark and luminous components of a 
RC, because very different mass models are thought to "produce" the same 
rotation curve 1 . The aim of this paper is to show that both these ideas are 
incorrect, i.e. that a) the characteristics of the observed RCs do not imply 
the existence of a fine-tuning in the halo and the disk structural parameters, 
but, on the contrary, these parameters take all the physically possible values 
and still produce smooth RCs, and b) "physically" different mass models yield 
detectably different RCs, provided that the observational data fulfill well de- 
fined high-quality requirements. 

We model the rotation curves of a galaxy from the distribution of the luminous 
and Dark Matter. We do not aim at the determination of the true Dark Matter 
density profile, issue that has to be dealt with by a different approach. Instead 
we assume the most popular density profiles, and we investigate the character- 
istics of the corresponding RCs. In detail, we study a pseudo-isothermal (PI) 
halo profile, suggested by observations (Persic et al. 1996), and a Navarro, 
Frenk & White (NFW) halo profile (Navarro, Frenk, & White, 1997) emerged 
in CDM N-body simulations. 

In section 2 we present the rotation curves derived for a PI halo profile (2.1), 
and for a NFW halo profile (2.2), and we discuss their characteristics. In sec- 
tion 3 we investigate the possibility of uniquely deriving, from an high-quality 
rotation curve, the actual mass model. Finally, in section 4 we summarize the 
results. 



Due to the results of this paper, we find no reason to explicitly identify the 
supporters of these ideas. 



2 



2 Generic mass models produce smooth, featureless Rotation Curves 



The total mass of a disk galaxy can be divided into two independent com- 
ponents, namely the luminous disk and the Dark Matter halo (neglecting the 
bulge is irrelevant in our study). These components contribute separately to 
the total rotation velocity: 

V 2 (r) = V d \r)+V h \r) (1) 



For the luminous matter, we consider the stellar disk, whose surface density 
profile is found as: 

E(r) = Z e- R / R ° (2) 



where S is the central surface density, and Re is the disk scale lenght 2 . The 
disk component of the velocity is (Freeman, 1970): 

Vlir) = u(r) (3) 



in which M D is the disk mass and u(f) = f 2 /2 • (IqK — IiKi)\ f / 2 , with I K 
and I\K\ being the Bessel functions of order and 1, evaluated at f/2. 
For the DM halo, we investigate both the PI profile and the NFW profile, 
(see figure (1) for comparison), with the halo velocity component being: 



VflO = ^ (4) 



where M H (f) is the dark halo mass inside the radius f. 

The contribution of the luminous matter to the total rotation velocity is deter- 
mined by two parameters: the disk scale lenght Rd, measured in each galaxy, 
and the total mass Mo, that can be obtained, at least in principle, from the 
disk luminosity, but since its determination is affected by large uncertainties, 
it is usually considered as an unknown quantity. The dark component of the 
rotation velocity has two free parameters, the central halo density and the 
halo characteristic scale lenght, that fully determine the DM distribution. 



2 It's worth noticing the observational relationship between the circular rotation 
velocity (V 3 = V(3R D )) and R D : R D ~5 (||) kpc (Persic et al. 1996). 
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Fig. 1. The normalized PI and NFW halo density profiles. 
2.1 Plhalos 



The DM density distribution of pseudo -isothermal (PI) halos is given by: 
Po 



p(r) = 



1 + {r/Rcf 



(5) 



For convenience we define: f = r/Rp as the galactocentric distance, and Rq = 
Rc/Rd as the core radius of the halo, both in units of the disk scale lenght. 
We take Rc as a free parameter and we consider its possible variation range 
between and infinity. Rc < 1 implies a mass distribution very similar to that 
of the NFW profile (see subsection 2.2). On the other hand, cases with Rc > 3 
are scarcely distinguishable from one another, since most of the observational 
data do not reach regions of disks beyond f = 3. 

The mass of the dark halo as a function of the normalized radius is obtained 



as: 



M H (f) = 4nGp Q R z c Rl 



Rc 



— tan 



Rc, 



(6) 



We define: 



£>3 



Rc 



— tan 




(7) 
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The total rotation velocity of the galaxy can therefore be written as: 



GM D G4np R 6 D 



V\f) = v{f) + " u " A(f) (8) 

riD tXD 



It is useful to perform the following transformation: 



_ GM D _ AnGpo g| _ ~ 3 



i . _i / i 



— tan" 



The parameters a and (5 are directly proportional to the fraction, at Rd, of 
the luminous and dark mass respectively. We define V± = V 2 (f = 1)). The 
radial profile of the normalized rotation velocity is then given by: 

V 2 (r) 

= a u(r) + [3 A(f) (9) 



a, f3 and A determine the shape of the curve, while V± determines its am- 
plitude. Let's stress that, as concerning the results of this section, the actual 
value of the constant V\ is irrelevant since, as a normalization parameter, it 
doesn't affect the RC shape. From (9) it follows that: 

az/i+/3Ai = l (10) 



where i>\ = v{r = 1). We notice that av\ and /3Ai represent the fractional 
contributions to the total (normalized) velocity at R D due to the luminous 
and Dark Matter respectively; for simplicity we define Jdm = /?Ai- 
We now rewrite the radial profile of the normalizated rotation velocity as 

772 ~ \ l ~ J DM ) H JDM ~: U-LJ 

V 1 Vi Ai 



under the following ranges of variations: 
< f DM < 1 < Rc < 3 

In this way, we split the original three free parameters into a 2+1 combination: 
2 determine the RC profile and 1 sets the RC amplitude. We allow the free 
parameters related with the profile to vary over a very wide range of values, 
and we plot the corresponding curves as functions of radius f. 
In figure (2) we present typical rotation curves we obtain for a variety of cases, 
including those with almost no Dark Matter in the central regions and those 
strongly halo-dominated, as well as those with large and small core radii. The 
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Fig. 2. PI rotation curves, for different /dm and normalized Rc- The curves rep- 
resent the total rotation velocity (thick line), the disk-component velocity (dotted 
lines) and the halo-component velocity (dashed lines). 
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Fig. 3. PI rotation curves with a second-order anomaly, indicating a sharp transition 
between a disk-dominated to a halo dominated region. 
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curves are plotted for < f < 4, covering the radial range of the observations. 
We realize that the variations of the parameters affect only the RCs' slopes: 
as foM increases, the curves become progressively more halo-dominated and 
therefore steeper, and the effect is enhanced for increasing R c . However, as the 
parameters vary, there are no peculiar features in the RCs, no flexes indicat- 
ing a discontinuity between a disk-dominated and a halo-dominated regime. 
In fact, the transition between these two regions appears smooth, and does 
not cause irregularities in the curves, at least in the zones actually mapped 
by data. No fine-tuning in the amount of the Dark Matter or in its density 
distribution parameters is required to obtain a smooth and featureless RC, 
like, for instance, those of the URC (Persic et al. 1996; see also Appendix A). 
Notice however that we can build a few peculiar models in which the transi- 
tion between a disk-dominated to a halo-dominated region makes itself visible 
in the RC profile, locally modifying the RC shape in a significative way. In 
figure (3) we show four of these cases. Anyway, they are very extreme, in that 
these galaxies have unrealistically low amounts of Dark Matter, Jom < 3%, 
and very large core radii, Rc > 6. Moreover, the flex in the curve, although 
quite significant, appears at galactocentric distances that are beyond the usual 
outermost radius reached by observational data. It is worth to point out that 
these peculiar RCs have the same profile often thought to be the standard 
outcome of the luminous and Dark Matter distributions in disk galaxies. In- 
stead, to be detected in the RCs, these features require a careful setting of the 
free parameters. The general properties of the dark and luminous matter that 
we observe in spirals make sure that, if a fine-tuning is needed to reproduce 
an observed RC, it is to feature irregularities, not smoothness. 



2.2 NFWhalos 



We now adopt the Dark Matter distribution of NFW halos (Navarro et al. 
1997): 



where the density is normalized to the critical density of the Universe, and 5c 
is the central density contrast (with respect to the external background). Rs 
is the halo characteristic lenght scale, and again we define: Rs = Rs/Rd an d 
f = r/R D . The mass of the dark halo inside a given radius is now: 



p(r) _ 5 C 



(12) 



Pent (r/Rs) (l + r/Rs) 2 



M H (r) = AT: Pcrit 8 c R%R% In 1 + ^ 



Rs 



r 



) 



(i + Rs/f) 



l 



(13) 
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After defining, in replacement of equation (7): 



r 



In 1 



RsJ (1 + Rs/f) 



(14) 



and using equations (1), (3) and (4), we write the expression for the total 
rotation velocity: 



2 GM D GAnp cr it5c R D 
V (r) = — — v(r) H - uiyr) 



R 



D 



R 



D 



(15) 



Similarly to the previous subsection, we set: 



a = 



GM D _ 4irGp crit 5 C R D _ 53 

P = ^ 779 > ^1 = K f 



R D V { 



2 ' 



RdV{ 



In 1 + ^~ 

. V Rs. 



l + Rs. 



and we get the normalized total velocity profile: 
V 2 (r) 



a u(r) + (3 uj(r) 



(16) 



with the constraint given by equation (10); finally, we obtain the parametrized 
expression for the rotation curve: 



v 1 — J DM) r J DM 



(17) 



Again, the two parameters f^M and Rs completely characterize the rotation 
curve profile. In figure (4) we show the NFW RCs for the same values of the 
parameters as the PI case, and we find a very similar behaviour: in general, 
NFW + disk curves are even smoother and flatter than the PI + disk ones; 
also in this case the emerging RCs have no second-order irregularities associ- 
ated to a luminous-to-dark regime transition, at least in the regions contained 
inside 4 disk scale lenghts. 



In the region contained inside the stellar disk edge, general PI and NFW 
mass models predict, for all values of the model free parameters, RCs that are 
smooth and featureless, i. e. with no flexes and inversion of the first derivative. 
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Fig. 4. NFW rotation curves, for different foM and normalized Rc- 
3 Disentangling a high-quality RC 

RCs are smooth and featureless when considered at small scales, and they 
show no individual charecteristics, apart from a different value of their "large- 
scale" slope. This is evident in the model RCs of the previous section and in 
the observational scenario provided by the URC. 

The further question is: how well does a rotation curve yield the true mass 
model of a galaxy? Moreover, is the information we obtain from the observa- 
tional data unique and model-independent? In other words, can we distinguish 
between different Dark Matter amounts and distributions by analyzing a ro- 
tation curve? In this section we will answer these questions for the PI case, 
which presents a number of difficulties, since it is partially degenerate with 
respect to its structural parameters f^M and Rq- The NFW case will not 
be considered here in that, having one free parameter less than the PI case 
and a known radial behaviour of the velocity profile, any result we find in our 
investigation for the PI will apply also to the CDM profile. 
Let's assume then that galaxies are composed by an exponential thin disk 
surrounded by a PI dark halo, and estimate the precision with which we can re- 
solve the luminous and the Dark Matter distribution from the rotation curve. 
Notice that the Universal Rotation Curve and the individual RCs are well 
fitted by a PI + disk model (Donato, Gentile & Salucci 2004; Gentile et al. 
2004). 

We build a family of 25 reference curves that simulate "observed" RCs of real 
galaxies, each one consisting of 25 data points, evenly spread between and 4 
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Fig. 5. The 25 reference "observed" rotation curves. 

disk scale lenghts, and is characterized by a set of two parameters in the en- 
semble: f DM = {0.1,0.3,0.5,0.7,0.9} ® R c = {0.1,1.0,2.0,3.0,4.0}. To each 
point of these curves we assign "observational" errors of the order to those of 
high-quality published RCs (Persic & Salucci 1995): ey = 0.02 is the averaged 
uncertainty in V(r), and = 0.05 in the RC slope. These curves are shown 
in figure (5). 

We investigate each of the "observed" RCs with a mass model, fitting the data 
to obtain their structural parameters. We remind that a successful strategy 
is accomplished by means of three suitable free parameters, namely the dark 
halo core radius Rc, the relative contribution of the Dark Matter to the to- 
tal velocity J'dm and the absolute value of the rotation velocity at a distance 
from the centre of one disk scale lenght V\. Of these parameters, two affect 
the shape of the RC, and one, V±, affects its amplitude. If the observed RC is 
of high quality, then it's possible to set the value of V\ directly from the data, 
reducing the number of free parameters. In this case, 

V 1 = V(l) ±C(10" 2 ) (18) 

Let us caution that this cannot always be done: we define a high-quality RC 
as a set of kinematical data for which such a parameter reduction is possible. 
Therefore, the reference "observed" RCs are fitted with the following mass 
model (see equation 11): 

VL d {f) = V\l) [(1 - f DM ) lL {¥) + f DM lDM {f] Rc)} (19) 

where j L (f) = v{f)jv\ is known, and j D M(r;Rc) = A(f; R c )/\i{Rc)- 
The measure of the distance between a model curve and an "observed" RC is 
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Fig. 6. Left: ^{V) for 3 reference curves, as a function of the model parame- 
ters: Rc € [0-2,4] (x-axis) and /da/ G [0.05,1] (y-axis). The "true" parameters 
are: (Rc = 1.0, /da/ = 0.3) (top panels), (Rc = 2, fDM = 0.5) (middle panels), 
(Rc = 3, /da/ = 0.7) (bottom panels). Right: Xtot f° r * ne same cases. 

given by a likelihood parameter, that we define as the sum of the x 2 computed 
on the velocity and on the RC slope: 

Xtot = X 2 [V] + X 2 [D] (20) 
where the terms on the right-hand side are: 




XD([fDM,Rc}mod) — 



J_ ( ^ n d Vmod(r n ; if DM, Rc)]mod) _ Jn_ dV obs (f n ) \ 

e 2 D —[\Vmod df V obs df J 

We denote with V mo d and V obs the velocities of the model curves and of the 
reference curves respectively. 
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The need of considering the contribution of the gradient of the velocity in the 
computation of the likelihood parameter is well explained by the examples 
shown in figure (6). In the left panels we show the value of the Xv om y> as 
a function of the values of the model parameters, for three different reference 
curves. The la contour of Xv includes also values of the parameters far from 
the "true" ones, unveiling a certain inability of the model to disentangle the 
mass components. With the information given by the gradient of the circular 
velocity and by the Xtot as defined in equation (20), the degeneracy of the 
different model curves largely disappears. In the right panels of figure (6) we 
show the resolution reached now: we find that the Xtot r i ses sharply to a 
value of ~ 10 2 in an interval of Af DM ~ 0.05 and of AR C < 0.25. These 
values set the resolution scale of the fit. 

For each one of the 25 "observed" reference curves, we perform the fit with 
the model given by equation (19), considering a very large range for the values 
of its parameters. We quantify the "success" of the fit by the computation of 
Xtot f° r each model curve considered. We now perform a general analysis 
of the behaviour of Xtot-> computed for 10.000 couples reference "observed" 
RC - model curve. It is worth to consider Xtot as a function of the following 
suitable "distance" in the paramater space: D = (A 2 (R C /A) + A 2 (/ DM )) 1 / 2 , 
with A{x) = x moc i — x b s (notice that, since the range of variation of Rq is four 
times that of J'dm, we normalize the former by a factor of 4 to make the two 
contributions to D comparable). The result is shown in figure (7a). The red 
circles represent cases in which the distance is dominated by the variation of 
the core radius ( A 2 (i?c/4) > A 2 (f DM )), the blue triangles cases in which the 
distance is mainly due to the variation of the amount of Dark Matter. The 
straight lines represent the limit value of Xtot corresponding to 3c and la. 
Most of the points lay over the 3a limit, indicating that the two curves of each 
couple are well resolved and therefore distinguishable. 

The vertical line defines a distance in the parameter space of D = 0.25. The 
points in the range D < 0.25 and laying below the 3a limit, although represent- 
ing cases in which the model fails in fitting uniquely a RC, can be considered 
irrelevant and will not be counted in the statystics: in fact, the parameters 
defining the curves are so close to each other, that even if a RC has multiple 
fits, they differ one from the other by a negligible amount. 
In the interval 1 — 3a we find about 0.9% of the total cases, and below the 
la limit about 0.23%. The cases below the la limit have to be considered as 
failures of the model in uniquely distinguishing the Dark Matter distribution, 
and those below the 3a limit as partial failures: therefore, studying our sim- 
ulated "observed" RCs we find that there is a total possibility of only about 
1% not to resolve (completely or partially) the mass model. 
Let's recall that we adopted the observational errors that are typical of about 
the top 50% of the observed RCs available today, for which, through the ap- 
plication of this procedure, there are very good possibilities to disentangle the 
mass distribution. However, it is possible to obtain an observed RC with an 
error even smaller than the one we adopted: for about the top 10% of the 
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Fig. 7. Xtot °f the 25 reference "observed" curves mapped with 400 model curves, 
as a function of the distance in the parameter space. Red points (circles) are for 
distances dominated by the variation of Rc, blue points (triangles) are for distances 
dominated by the variation of /dm- The straight lines represent the la and 3a limits. 
Panel a) e v = 0.02 and e D = 0.05. Panel b) e = 0.01 and e D = 0.03. 
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RCs available today we can consider observational errors of ey = 0.01 and 
6d — 0.03, and the number of these curves is going to increase in the future. 
The same analysis performed with these errors yields a Xtot vs D relation as 
shown in figure (7b). 

Now we have only about the 0.15% of cases in the interval 1 — 3a, and only 
the 0.02% is below la; the worst case among these sets the lower limit for 
the resolution of our model: (A (/dm)) max = 0.25 and (ARc)max = 1-2Rd, 
which is still a reasonable uncertainty. 

Therefore, considered the very small amount of cases below the 3a limit, 
we can state that the average resolution that we reach corresponds to the 
step with which we select the model parameters, i.e. (A(/om)) = 0.05 and 
(AR C ) = 0.2R D . 



4 Summary and conclusions 



We studied a model galaxy composed of a stellar disk embedded in a dark halo, 
and we constructed a family of rotation curves by varying the two parameters 
that affect the RC shape, i.e. the fractional contribution of the Dark Matter 
to the normalized rotation velocity at a reference radius, and the dark halo 
characteristic scale lenght. We investigated two different profiles for the distri- 
bution of the Dark Matter, i. e. a pseudo-isothermal and a Navarro, Frenk & 
White profile. In both cases we find that no fine-tuning is needed to reproduce 
the characteristics of the observed galaxy rotation curves, i.e. the evidence 
that they are smooth, with no irregularities like flexes in correspondence of 
the transition from the disk-dominated into the halo-dominated region. 
Eventually, irregularities do show up in the PI case, but only for extreme, 
tuned values of the Dark Matter content and distribution (foM < 0.1 and 
Rc > 6), and only beyond the disk edge. Since a similar result holds also for 
the NFW profile, we can conclude that, not only no fine-tuning is needed to 
reproduce the smoothness of the rotation curves, but the case is the opposite: 
a fine-tuning is required to obtain a "marked" curve. 

Secondly, with the purpose of investigating how effective is the process of 
mass modelling the RC, i.e. what is the precision with which we can infer a 
Dark Matter profile from the kinematics, we performed a statistical analysis 
on a large number of curves, assuming that the Dark Matter distribution is 
well described by a PI halo. We studied 25 reference rotation curves, repre- 
senting actual observed RCs, by means of a suitable mass model. We applied 
a maximum-likelihood method to assign a fitting model curve to each "ob- 
served" RC, and we found that in the 99% of cases we resolve the RC and 
distinguish between different Dark Matter distributions, while in the 1% of 
cases we found a degeneracy in the parameters. For this reason we can state 
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that, disposing of observations having a good accuracy and data for a suffi- 
ciently high number of radii, the distribution of the Dark Matter inside an 
observed galaxy can be determined from the RC. 

Of course, these results refer to the adopted profile used to describe the dark 
halo, for which the objection of the non-uniqueness of the mass modelling 
was raised. Also, our results don't minimize at all the importance of comple- 
mentary methods for determining the mass components of galaxies, like those 
coming from lensing, or those used to measure the disk masses (see Verheijen 
2004). 

With respect to the more general issue of determining the Dark Matter distri- 
bution in galaxies, several different profiles have been proposed to shape the 
dark halo, and the debate about which one is the favoured is still open. The 
coming of new-generation telescopes, that will provide high-precision observa- 
tions of the galactic dynamics, together with an improvement in the resolution 
of numerical simulations, will surely shed light on this topic. However, a deeper 
theoretical understanding of the structure of Dark Matter halos may still be 
required. 



A Universal Rotation Curve 
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Fig. A.l. The normalized Universal Rotation Curve. See Persic et al.1996 for details. 
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